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o 

CN One of the oldest ideas for unifying gravity and electroniagnetism is to con- 

bX) sider general relativity in five dimensions with one dimension curled up into 

a small circle. In other words, one studies vacuum solutions in five dimen- 
sions which asymptotically approach M^ x S^ at infinity. This is known as 



'^ Kaluza-Klein theory jH [2]. We have already discussed the simplest black 

hole solution in this theory: the product of Schwarzschild and a circle. In 



this chapter we discuss more general Kaluza-Klein black holes. We begin by 
Q^ describing a surprising property of the total energy in this theory. We then 

^ discuss black holes which are invariant under translations around the circle, 

and then drop this restriction. 



J^ 1 Energy in Kaluza-Klein theory 

1*:^ The total energy of an asymptotically flat five dimensional spacetime is given 

l/~j by a direct generalization of the standard ADM formula in four dimensions. 

(--.! Consider the metric on an asymptotically flat spacelike surface and let x* 

^D be asymptotically Euclidean coordinates. Then 

Tl gij = 6ij + hij (1.1) 

> 

r""! and the total mass is given by: 

"2 1 / ■ 

a M = ^^-^ j {hijj - hjj,,)dS' (1.2) 

where the integral is over a surface at large r. In Kaluza-Klein theory, an 
identical formula holds, but the integral is now over a surface with topology 
S"^ X S^, where the S'^ has radius r and the S^ has length L. The Kaluza- 
Klein vacuum, the product of a circle and four dimensional Minkowski space, 
has M = 0. 
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With standard asymptotically flat boundary conditions, there is a positive 
energy theorem which states, roughly speaking, that the solution with lowest 
total energy is Minkowski spacetime. More precisely, recall that the domi- 
nant energy condition states that T^j^t^'tg ^ for any two future directed 
timelike vectors t^,t2- One can prove the following |31S]: 

Positive energy theorem: Consider any nonsingular, asymptotically 
flat initial data set satisfying the dominant energy condition. Its total energy 
cannot be negative and vanishes only for flat Minkowski spacetime. 

The condition that the initial data be nonsingular is required to rule 
out, e.g., the negative mass Schwarzschild solution. This result holds in all 
dimensions I? > 4 [5J , but only if the spacetimes asymptotically approach 
M with the flat Minkowski metric. 

If a five dimensional spacetime asymptotically approaches M^ x S^ at infin- 
ity as required in Kaluza-Klein theory, then this theorem does not hold [B]. 
To see this, consider the five dimensional generalization of the Schwarzschild 
metric: 



ds'^ = - (l - ^] dt"^ + fl - ^\ dr^ + r^{d9^ + sin^ 9dn'^) . (1.3) 
We now analytically continue t ^ ix and — t- tt/2 + it. The result is 



2\ -1 

„2 _ / 1 '0 \ ,2 , / 1 ''O \ J 2 , 2, 



ds' = {l-'-^\ dx'+[l-^] dr'+ r'i-dr' + cosh^ rdn') . (1.4) 



Since (1.3) is Ricci flat, so is (1.4). To investigate the behavior near r = rg, 
set r = ro + p^/2ro. Then near /o = 0, the (x, r) part of the metric becomes 

dp' + p'^. (1.5) 

So if X is periodic with period 27rro, then p = (r = ro) is just a smooth 
origin of a rotational symmetry. Now consider the metric on the r = 
surface. This spatial metric is asymptotically flat, and since x is periodic, 
it satisfies Kaluza-Klein boundary conditions. What is the total energy of 



this solution? From (1.2) the energy comes from the 1/r correction to the 
metric, but since there is no 1/r correction, the total energy vanishes! This 
shows that M4 x S^ is not the only zero energy solution in the theory. 



What is the interpretation of the solution (1.4)? Since the x circle smoothly 
caps off at r = ro, there is no spacetime for r < ro- This is called a "bubble of 
nothing" , in analogy to a bubble of true vacuum which can nucleate inside a 
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false vacuum in field theory. This "bubble of nothing" expands rapidly until 
it is moving close to the speed of light and eventually hits null infinity. 

Not only are there nontrivial solutions with zero energy, one can find 
vacuum solutions with arbitrarily negative energy ^ . Consider time sym- 
metric initial data, i.e., initial data with zero extrinsic curvature just like the 
r = surface above. The only constraint on the spatial metric is that the 
scalar curvature vanishes. We now use the fact that the Reissner-Nordstrom 
metric has zero scalar curvature since the stress tensor for a Maxwell field 
in four dimensions is traceless. Analytically continuing t — )• ix and q ^ iq 
in the Reissner-Nordstrom metric yields 

ds2 = U{r)dx^ + U-^{r)dr^ + r^dQ'^ (1.6) 

with 

U{r) vanishes when r = r_|_ = ttt, -|- ^Jm? + q^, so we restrict the r coordinate 
to r > r_|_ . To avoid a conical singularity at r = r+ , we identify x with period 

471- 27rri 

U [rj^j r^ — m 

(Note that we are interested only in the metric and not the Maxwell field that 
is also a part of the usual Reissner-Nordstrom solution. We are constructing 
initial data for a five dimensional vacuum solution. In particular, q should 
not be thought of as an electromagnetic charge.) 



One can easily evaluate the energy (1.2) for this initial data and find 
G^M = mL/2. The parameter m can be positive or negative in this con- 
struction since r^ > and the initial data remains nonsingular. Thus we 
have negative energy solutions. (The special case tti = is precisely the ini- 
tial data for the "bubble of nothing" solution discussed above.) By taking 
m — )• — oo and q ~ (— m)^'^, one can make M arbitrarily negative keeping 
L constant. 

Given that there are solutions with arbitrarily negative energy in Kaluza- 
Klein theory, one might wonder why it is still taken seriously as a physical 
theory. The main reason is that all of the solutions we have discussed which 
violate the positive energy theorem have the property that spinors must be 
antiperiodic around the circle at infinity. This is a problem for the following 
reason: One can expand any five dimensional fermion in a Fourier series 
around the S^. The result is an infinite tower of four dimensional fields 
with masses which are multiples of 1/L. Since we want L to be small, the 
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only light fermion is the zero mode. In a space in which all femions are 
antiperiodic, there is no such zero mode. Hence all fermions will be very 
massive which contradicts the fact that we observe light fermions in nature. 
A realistic theory requires that we supplement the boundary condition with 
the requirement that fermions are periodic around the circle at infinity. 
With this added condition, one can prove a positive energy theorem [5]. All 
of the solutions we discuss in the remainder of this chapter will satisfy this 
condition. 



2 Homogeneous black hole solutions 

If a five dimensional spacetime is invariant under translations around a small 
circle, it can be viewed as an effective four dimensional spacetime coupled to 
certain matter fields. We begin by constructing the four dimensional theory 
which governs these fields. Letting y be a periodic coordinate around the 
circle with period L, we can write the five dimensional metric in the form 

ds^ = e-^'^/^{dy + 2A^dx^'f + e^'t'^^g^^dx^dx" (2.1) 

where (p,A^,g^iy depend on x^ but not y. The reason for the strange factor 
of \/3 in the definition of the scalar will become clear shortly. 
Note that under a coordinate transformation, y = y + 2\{x^), 

dy + 2A^dx'' = dy + 2A^dx^' (2.2) 

where ^^ = ^^+9^A. So A^ is like a four dimensional Maxwell potential and 
part of the usual five dimensional coordinate transformations are ordinary 
gauge transformations for A^. 



For the metric (2.1), y— ^ = e^'^l^'^yj^^ (independent of A^. If we 
evaluate the Einstein action for metrics of this form, we find (dropping a 
surface term): 

^=T7rV [dyd^xy^'R 
levrGs J 

g \R-2{V(l)f-e-^^^F'^] (2.3) 



1 

WttGa 



d'^x 



where F^i, = 2V[^A^] is the Maxwell field associated with A^, and G4 = 
G^/L. 

Thus, for metrics which are independent of y, the five dimensional Einstein 
action reduces to the action for four dimensional general relativity, a Maxwell 



field, and an extra scalar field (j). The strange factors in (2.1 ) were chosen so 
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that Qfj^ij has the standard four dimensional Einstein action. The equations 
of motion which follow from this action are 



y^/g-2V3<^^/.^\ Q (2.4) 

VV+— e-2^^F2 = (2.5) 



2g-2V3<^ p^^pp _ 1 g^^ ^-2V3<Pp2 (2.6) 



The five dimensional mass, computed from (1.2), agrees with the standard 
mass computed from the four dimensional metric g^,^. 

One simple class of solutions to these equations is </> = constant, F^jy = 0, 
and g^u is any four dimensional vacuum solution of Einstein's equation. This 
is just saying that if a four dimensional metric, dsl, is Ricci flat, then its 
product with 5*^ 

dsl = dsl + dy^ (2.7) 

is also Ricci flat. In particular, if dsl i^ ^^^ Schwarzschild metric, we get the 
black string discussed in the previous two chapters. 

The above relation between the four and five dimensional Newton's con- 
stant helps explain a puzzle about black hole entropy. In five dimensions, the 
entropy is proportional to the three dimensional "area" of the event horizon. 
Since this clearly increases with the length of the circle, one might think that 
the five dimensional black string can have many more microstates than the 
four dimensional black hole. But in fact the five dimensional entropy is the 
same as the four dimensional entropy: 

A5 A4L Ai 



4G5 4G5 4G4 



^4 • (2.8) 



2.1 Nonrotating charged black holes 

We have seen that neutral black holes can be obtained by simply taking a 
product of the Schwarzschild solution and a circle. The same is not true for 
charged black holes. These are not simply related to the Reissner-Nordstrom 
solution, since if F^^ ^ 0, it acts as a source for (p and </> cannot be con- 
stant. At first sight it seems difficult to find charged black hole solutions in 
this theory since the field equations involve exponentials of (p. However one 
can use a trick. From the five dimensional standpoint, one can generate a 
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Maxwell field by a simple boost in the y direction. Indeed, the four dimen- 
sional charge is simply the five dimensional momentum. So we start with 
the product of Schwarzschild and a line: 

ds^ = -(l-'^-^\ dp + (l - ^) "' dr^ + r^dn' + df . (2.9) 

We now boost this solution in the y direction 

t = t cosh a — y sinh a 

y = y cosha — t sinha . (2-10) 

Theifl 

, 9 / 2m cosh^ a\ ,9 / 2m sinh^ a\ ,9 
ds^ = -{1 dr + (1 + dy^ 



4?n, cosh a sinh a , , / 2?ti\ , 9 9 ,^9 
dtdy+ll J dr^ + r^dn^. (2.11) 

We now compactify y and put this into standard Kaluza-Klein form by 
defining 

e-^^/v^ = 1 + ?!!i!i!^ (2.12) 

r 

771 cosh a sinh a 

^* = - ^o ■ u2 ■ (2-13) 

r + 2m smh a 

In defining the four dimensional metric g^u one must remember to subtract 
the A^ term from gtt- Setting 

g = 2mcosh^a (2-14) 

the result is 

dr"^ 
g^udx^'dx" = -fdP + -^ + R^dn^ (2.15) 

where 

■^(^) = r 2^r~^r^ 11/2 ' ^'(^) = ^[^' + (*? - 2"^)^^]'^' • (2.16) 

[r^ + [q — 2m,)r\^i'^ 

This is a black hole with event horizon at r = 2m and singularity at r = 0. 
Note that the boost and dimensional reduction has not changed the location 
of the horizon or singularity. Unlike Reissner-Nordstrom, these black holes 
do not have a regular inner horizon. The singularity is spacelike, like the 
Schwarzschild solution. 

^ These boosted black strings are also subject to Grcgory-Laflamme instabilities [9]. 
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The total mass can be obtained by applying the general formula (1.2) to 



the five dimensional metric (2.11), but an easier approach is to compare the 
four dimensional metric (2.15) at large r with the Schwarzschild solution. 
Expanding the metric functions for large r, we have E? = r^ + r((/ — 2m)/2 + 
• • • and / = 1 — (g + 2rn)/2r + • • • . To compare with Schwarzschild, we take 
R to be our radial coordinate, so / = 1 — (5 + 2m) /2R = 1 — 2G4M/R. The 
total charge can be obtained by writing At = —Q/R + • • • . The results are 

G4M=i^, Q2^£(i_M_ (2.17) 

The extremal limit corresponds to the maximum possible charge for a 
given mass. This limit corresponds to m — )• 0, so one obtains 

Q ^ I , G4M ^ I ^ Q = 2G4M . (2.18) 

Recall that the extremal Reissner-Nordstrom black hole has Q = G^M so 
Kaluza-Klein black holes can carry twice as much charge. Actually, a Kaluza- 
Klein black hole cannot quite carry twice as much charge since the extremal 
limit is singular. The event horizon is at r = 2m,, so in the extremal limit 
{m = 0) the horizon becomes singular. Since the nonextremal black hole 
has a spacelike singularity, one might expect that in the extremal limit 
the singularity is either spacelike or null. However the following calculation 
shows that the singularity is actually timelike. Consider a radial null geodesic 
that hits the singularity. It satisfies 

dr 1 dr , , 

dt=j^q-^^^ (2.19) 

where the last expression is valid for small r. Since the ingoing radial null 
geodesic can reach the singularity at r = in finite t, it can lie entirely to 
the past of an outgoing radial null geodesic. This shows the singularity is 
timelike. If the singularity were null, then t would diverge as r — )• 0, just like 
geodesies approaching a horizon. 

What is the five dimensional description of an extreme black hole? Writing 



( 2.11 ) in terms oi q = 2m cosh a and taking the limit rn, — t- keeping q fixed 
yields 

ds'^ = -dt^ + dy^ + dr^ + r^dn'^ + - {dt - dyf . (2.20) 



r 



This is a five dimensional generalization of a plane fronted gravitational 
wave. Setting v = t + y,u = t — y, plane fronted waves take the form 

ds'^ = -dudv + dx'dxi + F{u, x') du^ . (2.21) 
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These metrics are Ricci flat if F satisfies a simple flat space Laplace equation: 

d'diF = (2.22) 

where i = 1,2,3 in five dimensions. The u dependence is arbitrary. The 
extremal black hole corresponds to a very simple plane fronted wave which 
is independent of u. It is a higher dimensional version of the Aichelberg-Sexl 
metric. 

We have discussed electrically-charged solutions. Magnetically-charged so- 
lutions can be obtained by a duality rotation: Setting 

*F^, = ^ e-^v^-^ e^/- Fp^ (2.23) 

the field equations are invariant under F — )• *F, (j) — )• —cp. (A similiar dual- 
ity rotation can be used to obtain magnetically charged Reissner-Nordstrom 
black holes.) So a magnetically-charged black hole has the same four di- 



mensional metric ( 2.15[2.16 ) but different matter fields. To emphasize that 
the solution is now magnetically charged we replace q with p, i.e., we set 
p = 2m cosh a. Then 



r 
A^ = P{l-cose) (2.25) 

where 

p2 ^ i^VF 1 ^ 2.26 

4 

P is the magnetic charge since Fg^j, = dgA^ = Psm0, so 

P = ^ f Fe^de d(j) . (2.27) 

There is no globally defined vector potential for a magnetic monopole, since 



if A were globally defined, Je2 dA = 0. The vector potential in (2.25) is not 
well behaved at 6 = it since J A^dcj) ^ even as the circle shrinks to zero 
size. This is sometimes called a "Dirac string". To avoid it, one can work 



with two patches on the sphere, using (2.25) on the northern hemisphere 
and setting A^f, = — -P(l + cosO) on the southern hemisphere. 

Let us reconstruct the five dimensional metric corresponding to the mag- 



2 Homogeneous black hole solutions 



netically charged black hole. Using the general formula (2.1) we get 

Note that the square roots that had been present in the four dimensional 
metric are gone. This spacetime still has a horizon at r = 2m and singularity 
at r = 0. We now take the extremal limit: rn, = 0. Remarkably, gu = —1, 
so the metric reduces to a simple product of time and a four dimensional 
positive definite metric. The four dimensional metric still looks singular at 
r = 0, but if we set x = 2//2-P and 



1 



2 pdr 



p = 2(pr)2 ^ dp'' = (2.29) 

r 

then near r = the spatial metric is 

2 

dp2 + ^ Udx + {I- cos e) d(l)f + de^ + sm'^ed(j)'^\ . (2.30) 

If X is periodic with period 47r, the quantity in brackets is the metric on 
S^ with radius two expressed as a Hopf fibration. So the metric is smooth 
at /O = (or r = 0) provided y has period SnP. This fixes the period of 
the Kaluza-Klein circle in terms of the magnetic charge P. This nonsingular 
solution is called the Kaluza-Klein monopole |10| \TT\ . 

It is a remarkable fact that Kaluza-Klein theory has a nonsingular mag- 
netic monopole solution since, e.g., Einstein-Maxwell theory does uoIq The 
spatial metric is called the Taub-NUT instanton and has topology M^. Asymp- 
totically, it looks like M^ x S^ and satisfies the usual Kaluza-Klein boundary 
conditions, but at small r the S^ combines with the 5^ of spherical symme- 
try to form an S^ which smoothly shrinks to zero at the origin of the four 
dimensional space. 

Since the four dimensional metric is the same as the electrically charged 
case, the extremal limit is singular. Thus we have an example of a four 
dimensional spacetime with curvature singularity that is resolved by lifting 
the solution to five dimensions. Conversely, one can see that the singularity 
in the four dimensional metric arises since gyy = at r = 0. In other words, 
the length of the circle that we are reducing along goes to zero there. 

We have seen that the extremal limit of both the electrically charged 

^ Of course the five dimensional metric is still a vacuum solution to Einstein's equation. The 
magnetic charge arises in the reduction to four dimensions. 
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and magnetically charged Kaluza-Klein black hole is singular in four dimen- 
sions. It turns out that black holes with both electric and magnetic charge 
(sometimes called dyonic black holes) have a nonsingular extremal limit 
with nonzero horizon area. We will not present it here, since the solution is 
considerably more complicated J12j and it is easily obtained from the more 
general rotating black holes discussed in the next section. 



2.2 Rotating Kaluza-Klein black holes 

To obtain a rotating, electrically charged Kaluza-Klein solution, one can 
repeat the construction in the previous section starting with the Kerr metric. 
In other words, one takes the product of Kerr and a line, boosts along the 
line and then compactifies the extra dimension. Rather than discuss this 
solution explicitly, we will jump ahead and give the most general known 
analytic family of black holes in Kaluza-Klein theory. These solutions are 
all stationary, axisymmetric, and invariant under translations in y. They 
depend on four parameters (m, g,p, a) which determine the mass M, electric 
and magnetic charges Q, P and angular momentum J. The five dimensional 
metric is: 



ds^ = ^{dy + 2A)2 -^{dt + B)2 + Hi(^ + dO^ + ^ sin^ Odcj? 



Hi' "^ ' ' H2' ' ' \ A ' ' Ae 

(2.31) 
where Hi, H2,Ag,A are all quadratic functions of r given explicitly by 



„2 , 2_2n , ^^^ o™^ , P (p-2m)(g-2m) 



Hi = r + a cos 9 + r{p — 2m) + 



p + q 2 



P 



2m{p + q) 
H2 = r^ + a^ cos^ 9 + r{q- 2m) + 



A/(g2^^4^^^2y(^^2T:4^ a cos 6* , (2.32) 



2,2 2/1 , „^„ o_^ , y {p-2m){q-2m) 



p + q 2 



+ — j^ — r A/f^^rri^^^^yTpr:^;;^ a cos i9 , (2.33) 

2m,[p + q) 

Ae = r^ - 2mr + a^ cos"^ 9 , (2.34) 

A = r^ - 2mr + a^ , (2.35) 
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and the l-forins A, B are given by: 

p — 2m\ Iq^ip^ — ^m?) 
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Q[r 



16m? (p 



-acos( 



H-^dt 



P{H2 + a^We)cose + 



p{q^ — ^m?) 



lQm?'{p + q)'^ 
X [{p + q){pr - m{p - 2m)) + q{p'^ - Am^)] asin^ 9] H~^d4{2.36) 



B 



(pq + 4m^)r — m,(p — 2m)(q — 2m) n 

pq z — -, r-^ a sin 



(2.37) 



2m{p + q)A0 

This comphcated solution was found by a solution generating technique 
which uses hidden symmetries of Einstein's equations with Killing fields 

After dimensional reduction, the five dimensional solution becomes a four 
dimensional black hole with the metric: 

A 



ds' 



.-^(.. + B)^ + yi?I^(^+..^ + -sm^ 



(2.38) 



The matter fields are the gauge field A given in (2.36), and the dilaton: 

H2 

Hi ■ 



-4<f>/V3 



2.39) 



The four parameters m,q,p,a appearing in the solution are related to the 
physical parameters M, Q, P, J through: 

p + q 



GaM 
GiJ 

p2 



'pq{pq + Am? 



q{q^ 



4{p + q) 



a 
m 



4(p + q) 

p{p^ — Ami?) 

4(p + q) 



(2.40) 
(2.41) 
(2.42) 
(2.43) 



The charge parameters Q, P were used already in writing the solution above. 
Note that q,p > 2m,, with equality corresponding to the absence of electric 
or magnetic charge, respectively. The solutions discussed in the previous 
section correspond to a = and either p = 2m or q = 2mi. 

These rotating black holes have inner and outer horizons at A = 0, i.e.: 



r± = mt 



± v^ 



m^ 



(2.44) 
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There are two qualitatively different types of extremal limits. Both yield 
black holes with smooth event horizons. The first is like Kerr: a = m. This 
is called the "fast rotation" case since G4J > QP. The other extremal limit 
is like the charged case discussed earlier. We take m — t- and a — t- 0, but keep 
the ratio a/m < 1 fixed. This is the "slow rotation" case since G4J < QP. 
In this limit, 

/ _i_ \ 2/3 
Q2/3 + p2/3^^i+P\ (2.45) 

so the mass of the black hole can be expressed in terms of the charges 

2G4M = (Q2/3 + p2/3)3/2 ^ ^2.46) 

The striking thing about this formula is that it is independent of the angular 
momentum J! As long as G4J < QP, one can add angular momentum to 
an extremal Kaluza-Klein black hole without changing the mass or making 
it nonextremal. Another surprising aspect about this extremal black hole is 
that the angular velocity of the horizon vanishes, even though the angular 
momentum is nonzero. This implies that there is no ergoregion. 

The black hole temperature vanishes in both extremal limits. The entropy 
takes a simple form. For the slowly rotating extremal solution {G4J < QP), 
the entropy is 



'5 = 2vr^/^^-J2 (2.47) 



'Q2p2 

while for the fast rotation case (G4J > QP), the extremal entropy is 



S = 27r J J2 - ^ . (2.48) 

The only difference is the overall sign inside the square root. When G4J = 
QP, the horizon area vanishes and the extremal solution is again singular. 
The solutions we discussed in the previous section have J = and only one 
nonzero charge, so their horizon area vanishes in the extremal limit as have 
already discussed. 



3 Inhomogeneous black hole solutions 

For the remainder of this chapter we shall focus on static vacuum five di- 
mensional solutions to Kaluza-Klein theory, but will drop the condition that 
the metric has an isometry associated to translations about the S^. With 
the isometry, the only static black hole solution which asymptotes to the 
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Kaluza-Klein vacuum is (2.28) J16tll7j. We further restrict our attention to 
solutions with no Kaluza-Klein magnetic charge, so that the only transla- 
tionally invariant static black hole solution is the black string. However, we 
shall see that the full space of static solutions without the circle isometry is 
very complicated indeedjj 

Without the S^ isometry, the Einstein equations lead to partial differential 
equations, and it is currently unclear how to find exact solutions. Pertur- 
bative techniques give a window onto certain solutions and we shall discuss 
these shortly. However generally we require numerical methods such as those 
discussed later in Chapter 10 to reveal the full and fascinating structure of 
the space of solutions. 

3.1 Localized black holes 

We will begin our exploration of the general static Kaluza-Klein black hole 
by taking a limit where we can neglect the fact that the extra dimension 
is compact. A very small black hole, where by 'small' we mean with radius 
much smaller than the size of the compact dimension L, will appear to a 
nearby observer to look approximately like a 5D Schwarzschild solution. The 
spherical symmetry of the black hole will be preserved near its horizon, but 
far away will be broken by the compactification of the spacetime. Conversely 
the asymptotic S^ translation invariance will be broken by the presence of 
this black hole which is localized at some particular position on the circle. We 
term such a solution a 'localized' black hole, and these were first discussed 
in detail by Myers [20] . 

Consider a 5D Schwarzschild solution. We may write this in isotropic 
coordinates as, 

dsk. = - (f^) ' dt' +(1+4!)' (dp' + ^'^^3) , (3.1) 

where the horizon radius is po, although we should take care that the coor- 
dinate position of the horizon, ph, is at ph = Po/2. Now of course very far 
from the horizon, for p^ p^ this solution becomes flat, 

dsl^f^ = {-dt^ + dr^ + r'^d^l + dy^) 

+ || {2df + dr^ + r'dnj + dy^) + 0(4) (3.2) 



3 



Much literature exists on this topic, and there are many details and interesting avenues we 
are not able to cover here. The interested reader is referred to two reviews on this topic 
I18II19| . although we note there has been significant numerical progress since these were 
written. 
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where we have chosen new coordinates, r and y to write this asymptotic 
behaviour, where, 

r = pcos6 , y = psinO , (3.3) 

so that /9^ = r'^ + y'^. The leading term looks just like the Kaluza-Klein 
vacuum, 

dsl^^ = -dt^ + dr"^ + r'^d^l + dy^ , (3.4) 

except that we cannot identify y ~ y + L as this is not compatible with the 



subleading term going as ~ 1/p in (3.2). 



Let us try to describe what happens to the metric away from the black 
hole more precisely. The method we use is based on the 'matched asymptotic 
expansion' approach of [211 I22j . We term the coordinates r, y the 'far field' 
chart. These are adapted to the asymptotic translation symmetry and y ~ 
y + L. The coordinates p,6 form the 'near field' chart and these are adapted 
to the approximate spherical symmetry near the horizon of a small black 
hole. These two charts must overlap. Let us take the extent of the near and 
far field such that, 

far field : r^ + y^ > (i _ a)^^ , near field : pl<p^ <{1 + A)e^ (3.5) 

where £ is a scale parametrically larger than po and smaller than L in the 
limit po/L —7- 0. We might, for example, take i'^ = poL. The constant A is 
0(1) and ensures the charts overlap so that < A < 1. 

Far from the horizon the black hole appears to be a localized source of 
mass, and being static, we expect it can be treated as a point mass source 
in Newtonian perturbation theory. For a point mass source at r = y = in 



the vacuum (3.4) one may solve the linear problem to obtain the metric in 
the far field as, 

ds^ ~ dsl^^ + $ (2 dt^ + dr'^ + r'^dnj + dy^) , (3.6) 

where the Newtonian potential is explicitly given as, 

$ = !^5^ ^.1^^^ (37) 

3Lr cosh 2ji:+ cos ^ ' ^ ' 



where M is the mass as computed in the earlier equation (1.2) and $ is a 
harmonic function on ( |3.4[ ) with a delta function source at the origin |20| . 
In the overlap region then p/L <^ 1 and we may expand this potential in 
the near field coordinates as. 
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and taking the leading two terms we write the metric in the overlap as, 

ds' c. dsl^ + ^ (1 + ^) (2 dt' + dp' + p'dnl) . (3.9) 

In the near field the matched asymptotic expansion method instructs that we 
consider perturbations of 5D Schwarzschild. However, for our discussion the 
perturbation we require can simply be thought of as a global scaling together 
with a scaling of time. Hence the near field solution remains Schwarzschild 
which we write as, 

+ (1 + %) (l + ^) (V + P^dnl) , (3.10) 

where q and Cy are constants, with |ct|,|cj^| <C 1, which give the required 
perturbative scalings. Then in the overlap, where p^ po the metric behaves 



similarly to the expansion (3.2) giving. 



2 ^, J„2 , Po /o j+2 , ,2 , 2 ,n2^ 



ds'^ ~ dsl^^ + ^ {2dt' + dp'' + p'dft: 



2p 



3) 



- ctdf + Cy {dp' + p'dnl) , (3.11) 

We see that in order to 'match' the far and near field asymptotic expansions 
we must identify, 

G,M=^-^pI, -lct = cy = '^. (3.12) 

We have matched the first two terms in these asymptotic expansions. To 



proceed further we see that the next term in (3.8) involves an angular de- 
pendence and cannot be matched simply by the Schwarzschild metric alone. 
One must perform static perturbation theory about the Schwarzschild met- 
ric, and match this in a multipole expansion. In addition, in the far field one 
must go beyond Newtonian order. This procedure is detailed in [21] which 
shows how to construct the metric in the near and far field order by order in 
the perturbation parameter po/L. The simple behaviour we have discussed 
above arises since the lowest multipole static perturbations of Schwarzschild 
(being spherically symmetric) can only be trivial global and time scalings 
due to Birkhoff 's theorem. 

Having related the far field parameter, the mass M, to the near field 
parameters po, ct and Cy, we may proceed to compute physical quantities. 
From the near field solution we may compute the horizon area. A, and surface 
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gravity, n, perturbatively in powers of p^/L, as, 

A ^ 2.^pI{1 + cyfl^ = 2.-pI (l + ^ + O (|| 

We know that the mass has leading behaviour M = ^-p^- However, we 
can compute the subleading correction to this using the first law at constant 
circle size L, so we have, dM = „3j dA. Then using the expressions for A 
and K above yields, 

and hence we see that for small localized black holes we have the behaviour, 



where the leading behaviour is that predicted by Myers PU], and the sub- 
leading corrections were computed in [ 22| [23] . 

Another interesting quantity to compute is the proper distance between 
the poles of the horizon along the axis of rotational symmetry, Laxis- One 
might imagine that as the black hole size is increased there is a corresponding 
linear response decreasing the length of the axis as Laxis — L {l — a^\ for 
some constant a > 0. This constant may be computed from our metrics 
above, and surprisingly one finds it precisely vanishes. There is no linear 
variation of Laxis, and this has been termed the 'Archimedes effect' |21j . To 
leading order, the geometry around the small black hole precisely expands 
to accommodate it. 

Moving beyond perturbation theory we expect that the localized black 
holes exist as solutions to the full Einstein equations at least for G^M <^ L^. 
Since the 5D Schwarzschild solution is dynamically stable we expect that 
the localized black hole solutions are similarly dynamically stable in this 
limit. An interesting question we shall return to later is what happens to 
this branch of solutions as one increases their mass for fixed L. Presumably 
when G5M ~ L^, so that their horizon size becomes of order the asymptotic 
circle size, they become strongly deformed away from the Schwarzschild 
geometry. One possibility is that the inter-polar distance Laxis actually goes 
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to zero at some finite mass, say M^, so that the branch of solutions ends. 
Another possibihty is that the Archimedes effect ensures that the geometry 
always accommodates an increasingly large horizon, and there is no upper 
mass limit to these solutions. We shall return to this question later, and 
find evidence that the former suggestion occurs. However, before we discuss 
this we must introduce a rather exotic static black hole, the inhomogeneous 
black string. 



3.2 Inhomogeneous black strings 

Earlier we discussed the simplest static black hole in Kaluza-Klein theory, 
the black string, which is a straightforward product of the Schwarzschild 
metric and a circle. Let us write this solution in Schwarzschild form as, 

„2 _ /i '^0\ ,.2 , 1 j„2 , „2 ,r.2 , j,.2 



ds^ = - M - ^ j dt^ + z—T^dr^ + r^dni + dy^ . (3.16) 

r 

We have seen that this solution is dynamically unstable to perturbations 
with wavenumber k on the circle less than a critical wavenumber kc — 
0.876/ro. Let us choose our extra dimension to have length L so that y ~ 
y+L, and take this length to precisely be the critical wavelength for marginal 
stability, L = 27r/kc. Being marginally stable this black string has an ex- 
actly static linear perturbation and Gregory and Laflamme realized that 
this might signal the existence of a new class of static solutions [24j . Gubser 
argued that this deformation does indeed lift to the full non-linear theory 
to generate an entirely new branch of solutions [25] . Since the perturbation 
explicitly breaks the translation invariance on the circle these static solu- 
tions have the horizon topology of the black string, but are inhomogeneous 
on the circle. We term them inhomogeneous black strings, in contrast to 



the solution (3.16) above which we now refer to as a homogeneous black 
string. Such inhomogeneous black strings were first discussed by Horowitz 
and Maeda as possible end states of the Gregory-Laflamme (GL) instability 
|26j . Interestingly, at least for 5D Kaluza-Klein theory it appears that they 
are all unstable, and therefore cannot serve as such an end state. 

Just as for the localized solutions, these inhomogeneous black strings may 
be constructed perturbatively. For the localized black holes the perturbative 
limit was when the horizon size was very small compared to L. For these 
inhomogeneous black strings the perturbative limit is when they are very 
weak deformations of the homogeneous marginally stable solution. We now 
give an overview of this perturbative construction, essentially following the 
approach of Gubser [25J . For this discussion we will perform an overall scaling 
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of the solution to fix rg = 1, so that L = 27r/0.876. We may then write the 
general metric with the isometries of the inhomogeneous black strings as, 

ds^ = -(l-^^ e'^dt^ + e^^ [y^J^^ + dy^\ + r^^'^'dnl , (3.17) 

where the functions A, B, C depend on r and also y due to the inhomogene- 
ity. We require that the horizon at r = 1 is regular and that asymptotically 
A,B,C^ which ensures the static Killing vector d/dt has unit normal- 
ization and the circle size is L. 

At linear order in perturbation theory we have the marginal GL mode. 
This takes the form, 

A = X a{r) cos kcH , B = Xb{r) cos key , C = X c{r) cos key , (3.18) 

where A is the perturbation parameter, and a, b, and c are functions of r that 
must be determined numerically by solving ordinary differential equations. 
Gubser's method gives a systematic way to compute the backreaction of this 
mode. Non-linear terms in the Einstein equations couple the various Fourier 
modes on the circle. The linear term squares to give a source at quadratic 
order, and since cos^ key = ^(1 -|- cos2kcy) one obtains contributions to the 
backreaction at O(A^) going as the constant mode and cos 2kcy on the circle. 
Similarly since cos^ kcV decomposes into cos key and cos 3kcy components, 
at cubic order O(A^) one has backreaction in these Fourier modes, and so 
on. Thus the full solution generated by the backreaction of the marginal GL 
mode takes the form, 

Aai^i cos kcV 
+X^a2,o +X^a2flCos2kcy 

A= +A^a3^i cos fccy -|-A^a3_3 cos 3A;cy (3.19) 

+A^a4^o +A^a4^2 cos 2kcy +A^a4^2 cos Akcy 

+ ... 

where an,m are functions depending on r, and B and C have expansions 
taking the same form, with coefficient functions hn^m and Cn,ni respectively. 
We denote the level (n, m) as being at order 0(A") and Fourier mode 
cosmkcy- The leading linear term is the level (1, 1) marginal GL mode, so 



that ai^i(r) = a{r) in equation (3.18) and similarly for bi^i and ci^i. At the 



level (n, m) one must solve the problem which we schematically represent 
as, C{an,m,bn,m,Cn,m) = <S where the left-hand side, £, is a homogeneous 
linear differential operator acting on the level {n,m) functions an,m, bn,m 
and Cn^rm and the right-hand side S is the source for this inhomogeneous 
linear problem. This source originates not only from the component of the 
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n-th power of the leading linear term ai^i, 61^1 and ci^i in the m-th Fourier 
mode, but also from various combinations of the intermediate backreaction 
orders. 

Solving these inhomogeneous linear systems is technical but straightfor- 
ward, and we will not give the details here. The level (n, m) linear system 
can be computed by numerically integrating ordinary differential equations 
provided all the previous orders 0(A), O(A^) . . . 0(A"~^) have been com- 
puted. In practice computing levels (1,1), (2,0), (0,2) and (3,1) is rather 
straightforward and that is all we shall require for our discussion. Numerical 
computation gives the following data. 



fee = 0.876; 



ai,i(l) = 


= ^1,1(1) = 


= -0.55, 


ci,i(l) = l 


a2,o(l) = 


= -0.28, 


^2,0(1) 


= 0.77, C2,o(l) 


«2,2(1) = 


= ^2,2(1) = 


= 0.34, 


C2,2(l) = -0.69 


«3,l(l) = 


= ^3,1(1) = 


= -0.24, 


C3,i(l) = 



0.80 



(3.20) 

which in principle specifies the solution to the levels (1,1), (2,0), (0,2) 
and (3, 1) if one integrates these data out from the horizon to infinity (al- 
though one would require more precision than the two significant figures 
given here) Ij For these levels the functions an,m , bn,m and Cn,m all expo- 
nentially decay asymptotically except for (2, 0) which has power law decay 
going as, 

, Soo „. 1 . Coo log r 1 

^2,0 = ^0{^), C2,0 = ^0{-), 

5oo = 0.41 , Coo = -0.12 . (3.21) 

The variation of the mass M, area A and surface gravity k for a non-uniform 
string takes the form, 

M = MGL{l + m2\^ + miX^ + ..) , A = Aql {l + a2X^ + a^X^^ + . . .) , 
k = kgl(1 + K2A2 + K4A^ + ...) , (3.22) 

where, 

0.876L2 ^ 0.8762L3 vr 

^- = ^;^' ^- = ^^' ^- = 0:876L' ^'-''^ 

and we have now rescaled our solution to have a circle size L again. Note 
that the marginal GL perturbation of level (1,1) does not contribute to these 
at order 0(A) as it decays exponentially away from the horizon, and being 
a harmonic perturbation on the circle, does not change the horizon area at 

■* Note that we have chosen the constant circle size 'scheme' of 1251 so that the wavenumber in 
(|3.19|l is unperturbed at higher orders in A. 
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this order. In fact we see the corrections to these quantities only arise from 
even powers in A. The quadratic variations are determined from the above 

m2 = 35oo - 2Coo = 1.45 

02 = 62,0(1) + 2c2,o(l) + \h^iilf + 6i,i(l)ci,i(l) + ci,i(l)2 = 2.90 

K2 = ai,i(l)- 61,1(1) = -1.04 (3.24) 

and thus in a similar manner to the small localized solutions, we have com- 
puted the properties of the inhomogeneous string in a perturbative limit. 

Let us consider the dynamical stability of the homogeneous black strings. 
Homogeneous strings with M < Mgl are unstable to GL perturbations with 
wavelength L. Indeed for small enough mass, M < McL/n for integer n the 
wavelengths L/2, L/3, . . . , L/n will also be unstable as the higher harmonics 
of the GL mode fit into the circle. For M > Mql the homogeneous black 
strings are stable as the unstable GL modes, having minimum wavelength 
(M/Mgl)L, cannot fit onto the circle. 

Now let us understand the stability of the weakly inhomogeneous solu- 
tions, namely those with A ^ 1, so that M is close to Mgl- We see that 
since 777-2 > these inhomogeneous solutions have mass greater than Mgl 
and hence coexist at the same mass as a stable homogeneous black string. Let 
us compare the areas of these two solutions. Consider moving an infinites- 
imal distance dX along the inhomogeneous black string branch starting at 
A = 0. Then the first law, dM = KdA/{%TTG^), for fixed circle size implies 
at order O(A^) that 27772 = ^2; which we see is consistent with the numer- 
ical data above. At order O(A^) it implies that 27n,4 — 04 — \n2a2 = 0. A 

homogeneous string has area Ah = ( — 'jj-^ I M^ for a mass M, and so we 
can compute the fractional area difference between the inhomogeneous and 
homogeneous strings for a mass M as, 

/S.A A — Ah 

—- = — - — - = (a2 - 27772) A^ + (04 - 2a2"72 + 8777^ - 27774) A"^ + O(A^) 

Ah Ah 

= -7772 (^2 + K2) A^ + 0(A6) ~ -0.59A^ (3.25) 

where the quadratic term vanishes by the first law, and the quartic term is 
precisely determined by the data we have given above. We arrive at the result 
that at fixed mass and fixed circle size, the weakly inhomogeneous string has 
a lower area than the stable homogeneous string of the same mass. Since 
area cannot decrease in a dynamical process this argument shows that the 
weakly inhomogeneous strings are dynamically unstable to perturbations 
that preserve their mass, and deform them into the homogeneous strings. It 
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has currently not been demonstrated whether this perturbative instabiUty 
can be seen in linear perturbation theory or only at higher orders. However, 
we believe it is likely to manifest as a linear dynamical instability, and we 
shall later assume this is the case. A reasonable supposition is that such a 
linear instability generates an evolution that ends at a stable homogeneous 
black string with similar mass. 

As with the localized solutions, an interesting question is what happens to 
the branch of inhomogeneous solutions as one deforms past the weakly in- 
homogeneous regime so that the perturbative approach above breaks down. 
In the next section we shall consider this. 



3.3 The space of static black hole solutions 

We have seen that the homogeneous black strings exist for all mass, with 
fixed circle size, but are stable only for sufficiently large mass, M > Mql ~ 
O^L'^/G^). At low masses, M <C L'^/G^, we have another branch of solu- 
tions, the localized black holes, which we expect to be dynamically stable. 
In the intermediate mass range for masses just above Mqj^ we also have 
the inhomogeneous black strings, which we have argued are unstable. We 
already see there is no uniqueness for static solutions to Kaluza-Klein theory 
for fixed mass and circle size. 

We have posed the question of what happens to localized black holes 
as they become large, and also what happens to the inhomogeneous black 
strings as they become increasingly inhomogeneous. Harmark and Obers 
|27j argued that the localized black holes might continuously connect to the 
inhomogeneous black string solutions proposed by Horowitz and Maeda. Kol 
used a simple Morse theory argument to deduce that the simplest scenario is 
for the localized black holes to continuously connect to the inhomogeneous 
black string branch found by Gubser [2H] . Furthermore he predicted that the 
two branches connect via a singular solution that mediates a topology change 
of the horizon. This is analogous to the topology change between hyperbolic 
and parabolic conic sections, which is mediated by the singular section taken 
through the apex of the cone. Indeed there is a static, spherically symmetric 
Ricci flat conical geometry, 

dsl^^ = da^ + -Q^ (d/32 - sin^ /3dt^) + -a^dQ^ , (3.26) 

which is singular at its apex a = 0. The base of the cone is the product of the 
2-sphere with the 2D de Sitter space (d/S^ — sin^ f3dt'^), with (3 e [0, 2tt] where 
/3 = and /3 = vr are Killing horizons with respect to d/dt. These horizons 
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are connected (although not smoothly) at the apex of the cone where the 
de Sitter factor is zero sized. Kol proposed that this gives a local model for 
the geometry of the topology changing solution at the singular point of the 
horizon. Moving away from the singular topology changing solution to the 
inhomogeneous black string branch would resolve the sphere so that it was 
finite in size at a = 0, in accord with the inhomogeneous strings having 
no exposed axis of symmetry, and the two components of the horizon with 
/3 = and /3 = vr would be smoothly connected. Conversely, moving to the 
localized black hole branch would resolve the de Sitter factor, so that the 
horizons at /? = and /3 = vr would not touch, and a = would be the 
exposed axis of symmetry. 

In order to make progress in understanding the extension of the localized 
black hole and inhomogeneous black string branches we must resort to nu- 
merical work. The inhomogeneous black strings were first computed in [2^ 
and the localized black holes in [301 [3T| l32] . Recently a standard approach 
to numerically construct static vacuum solutions has been developed [33] 
and will be discussed in Chapter 10. Unlike the previous methods employed, 
for example [29j , this approach is geometrically elegant being covariant and 
may be applied to problems depending non-trivially on as many coordinates 
as one likes. In practice these new methods work very well. The results pre- 
sented here for both inhomogeneous black strings and localized black holes 
are those computed in [33j . 

As we shall now see, the results are compatible with Kol's prediction of 
a topology changing merger of the localized black holes and inhomogeneous 
black strings discussed above. The numerical calculations have not yet been 
adapted to explore the singular potential merger point, and inevitably break 
down for both the localized and inhomogeneous string branches as this point 
is approached. We note that it is still early days for these methods, and we 
fully expect that future work in coming years will significantly improve on 
the solutions reproduced here. 

Let us begin in Figs. [L] and [2[] by showing the curves of area and sur- 
face gravity, plotted as the inverse temperature 1/T = 27r/K, against mass 
for fixed circle size L. Quantities are made dimensionless with appropriate 
powers of the circle size L. We find a very interesting behaviour for the 
localized black holes, namely that there is a maximum mass solution with 
G^Mmax ~ 0.17 L^. We note that since dM = KdA/SirG^ for fixed circle size 
L, then any extremum of M will similarly be an extremum of A. Hence the 
area against mass curve has a cuspy profile at M = Mmax- We also see from 
Fig. ^ that there is a minimum surface gravity (or equivalently tempera- 
ture) solution, at M = M^. The actual value of M^ is numerically close to 



3 Inhomogeneous black hole solutions 



23 



AIU 



1.2 



1.0 



0.8 



0.6 



0.4 



0.2 



Homogeneous Strings 



Localized Black 
Holes 




Localized Black Holes 



Merger point? 



Inhomogeneous Strings 



GL point 



100 



0.05 



0.10 



0.15 



GsM/L^ 



Figure 1. Plot of area A against mass M for fixed Kaluza-Klein circle 
size L for the homogeneous (dashed) and inhomogeneous (solid thin line) 
black strings and the localized black holes (solid thick line). The curves 
for the latter two are the numerical solutions of [33] ■ The inhomogeneous 
black strings and localized black holes are compatible with a merger at 
GsM^/L^ ~ 0.12, and there is a maximum mass localized solution, with 
G^Mmax/L'^ — 0.17. The g ray d otted line is the small localized black hole 
approximation in equation (3.15 1 and we sec the approximation is excellent 



for increasing mass up to M ^ M, 



that of Mmax, as is clear visually, but we emphasize that since the curve 
of 1/k, against M is smooth for the localized solutions about the maximum 
mass, as we expect and indeed see numerically, the point of maximum mass 
and minimum surface gravity cannot coincide, and M^ < Mmax- This leads 
to the peculiar conclusion that the localized solutions in the narrow range 
Mk < M < Mmax actually have positive specific heat. 

We confirm that the perturbative prediction computed above in equation 



(3.15 ) provides a good description of the localized behaviour for small mass. 



In both figures this approximation is plotted, and we see surprisingly good 
agreement for increasing mass up to Mql and even past this point to around 
M ~ Mmax- These approximations have no maximum mass, and hence 
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Figure 2. Plot of inverse temperature 27r/«; against mass AI for fixed circle 
size L for the same solutions as in the previous figure. The labels 'A' to 'I' 
denote solutions whose horizon embeddings are displayed in the following 
figure. Again we see consistency with a merger. We also note that there is a 
minimum surface gravity solution at M^j (near the label 'C') which is rather 
close to, but sl ightly less than Mmax- The perturbative approximation in 
equation (3.15) is also plotted and again very good agreement is seen for 



localized solutions with increasing mass up to around M ~ M„ 



GsM/L^ 



cannot agree past this point in the localized branch, and we expect non- 
perturbative effects to become important there. 

The inhomogeneous black string branch departs from the homogeneous 
one at the GL point M = Mgl- Moving away from this point the deforma- 
tion of the horizon becomes greater. The degree of inhomogeneity, measured 
by the ratio of the maximum and minimum 2-sphere radii of the horizon, 
increases monotonically. The most striking feature of these plots is that 
moving as far along the localized and inhomogeneous branches as possible, 
both branches appear to be compatible with joining at some mass M^ with 
G^M^, ~ 0.12L^. This is an indication that indeed these branches both end 
at the same singular topology changing solution. Other physical quantities 
that we have not displayed here also support this conclusion. 
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Figure 3. Figure showing the spatial geometry of the horizon for a number 
of locaHzed black holes and inhomogeneous black strings all with the same 
asymptotic circle size. The embeddings are labelled 'A' to T and correspond 
to the solutions annotated in the previous figure. We emphasize that these 
geometries are those actually found in the numerical solutions of [33j . For 
the localized black holes the proper length of the axis of symmetry is also 
depicted. 



In Fig. [S] we show several embeddings of the horizon for various locahzed 
and inhomogeneous string solutions. We see graphical confirmation that the 
most inhomogeneous strings found appear to be close to having the minimum 
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2-sphere of the horizon go to zero size. The locahzed black holes also appear 
to be close to the poles meeting as the axis of symmetry shrinks to zero 
length. Near the potential merger, the geometries appear to be similar away 
from the singular region. 

We note that by uncompactifying the solutions above, and recompactify- 
ing with circle size an integer multiple of L one can construct multi-centered 
localized black hole solutions although all such solutions presumably would 
be unstable. 

The metric for the inhomogeneous black strings has been tested to see 
explicitly that the cone geometry advocated by Kol emerges for solutions 
with M ~ M^ in the region surrounding the minimum radius of the hori- 
zon. Strings in various dimensions were found, and the appropriate cone 
geometries were indeed seen to emerge [Si]. Similar tests have not yet been 
performed on the localized solutions, but we believe it is very plausible that 
this cone geometry will again be seen. 



3.4 Stability of the static black hole solutions 

Having seen an elegant proposal for the space of static solutions and nu- 
merical evidence supporting it we now discuss their likely dynamical sta- 
bility. From the earlier Fig. [L] we see that for fixed circle size the local- 
ized solutions dominate the area for a given mass when M < Merit where 
G^Mcrit — 0.12L^, and for M > Merit the homogeneous black strings domi- 
nate the area. Thus at fixed mass (the 'microcanonical ensemble') we expect 
that the globally stable solution is a localized black hole for M < Merit and 
a homogeneous string for M > Merit- 

This does not tell us about local dynamical stability of the solutions. For 
example, the homogeneous strings with Mql < M < Merit are dynamically 
stable to small perturbations even though they do not globally dominate the 
area, and can be deformed to a more stable solution with the same mass for 
a sufficiently large deformation. 

Let us now consider this local stability. We have argued that the small 
localized solutions are perturbatively stable and the weakly inhomogeneous 
black strings are perturbatively unstable. Furthermore it appears that these 
two branches should be thought of as one, joined together by a singular 
topology changing solution. If these assertions are correct the simplest sit- 
uation is that either i) the entire inhomogeneous branch and some of the 
localized branch are unstable to a single perturbation mode, or alternatively 
ii) the instability ends at some point along the inhomogeneous branch and 
so there are strongly deformed stable inhomogeneous solutions and all the 
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localized black holes are stable. For either of these scenarios to be true the 
instability should smoothly turn off as one moves in the space of solutions, 
and hence there will be a localized solution (for scenario i) or an inhomoge- 
neous string (scenario ii) with a marginal static perturbation. 

The localized branch has precisely such a candidate marginal static per- 
turbation as we shall now argue, indicating that scenario i) is realized. Since 
we are considering the microcanonical ensemble we require a static pertur- 
bation that leaves the mass invariant. (A static perturbation that changes 
the mass represents a change in the black hole, not a perturbation of a given 
black hole.) Take the perturbation of the metric for a localized solution to 
be that generated by moving infinitessimally along the branch of solutions 
(i.e. the perturbation is the tangent to the space of solutions there). For a 
generic solution such a perturbation will change the mass to linear order. 
However precisely at the point M = M^ax the perturbation will leave the 
mass invariant to linear order. This then gives the required candidate static 
perturbation that is both regular at the horizon and leaves the mass invari- 
ant jj Physically we may gain some intuition for this change in stability at 
the maximum mass solution by considering taking a small (presumably sta- 
ble) localized black hole and gradually dropping matter into it. Initially one 
may expect the mass will grow as it absorbs the matter. However, when it 
reaches the localized maximum mass solution with M = Mmax (and corre- 
spondingly maximal area) it can no longer remain a stable solution. Adding 
more matter will force the area to increase during the ensuing dynamics 
implying the solution cannot settle back to any solution near Af = M^ax 
since these all have lower area. 

Let us summarize our stability discussion. The numerical solutions imply 
that localized black holes with M < Merit are globally dynamically stable, 
while homogeneous strings with M > Adcrit are globally stable. We have 
argued that the simplest picture of linear stability is that the localized black 
holes are linearly stable for small mass, and remain stable as their mass is 
increased up to the value M = Mmax- Moving further along this localized 
branch these solutions become unstable to a single perturbation mode, and 
this instability continues through the topology change, and all the way along 
the inhomogeneous black string branch. A further comment is that whilst 
dynamically unstable solutions may be interesting in terms of understand- 
ing the moduli space of static vacuum solutions, presumably they play no 
physical role. Whilst we have had to resort to sophisticated numerical work 
to elucidate the full structure of solutions, from a physical perspective the 

^ We thank Harvey Reall for an important discussion clarifying this issue. 
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interesting (i.e. stable) solutions are the simplest ones. The homogeneous 
black strings we know analytically, and the stable localized solutions with 
M < Mmax appear to be rather well described by the perturbative construc- 
tion we gave earlier. 

Sorkin has shown this picture changes remarkably in higher dimensions 
|35j . Treating the dimension as a continuous parameter D, again with a sin- 
gle compact Kaluza-Klein circle, so that previously in this chapter we have 



considered only D = 5, then one finds that the sign of m2 in equation (3.22), 



and correspondingly that of the area difference in equation (3.25), changes 
above a critical dimension D^, ~ 13.5. Now weakly inhomogeneous strings 
have masses M < Mgl, and are expected to be stable and the end state 
of the GL instability of homogeneous strings with mass just below Mgl- 
Assuming the localized solutions and inhomogeneous strings still merge at 
some mass M^, then the simplest picture is that for M < M^ the localized 
black holes exist and dominate the area for fixed mass over the homogeneous 
strings, and for M* < M < Mql the inhomogeneous strings exist and dom- 
inate the homogeneous ones. An interesting question is whether there are 
inhomogeneous black strings that are stable for D < D^. We have argued 
for D = 5 they are likely to all be unstable. However, for D close to, but less 
than D^, continuity implies the strongly inhomogeneous black strings will 
remain stable, even though the weakly inhomogeneous ones will become un- 
stable. Whether stability persists for strongly inhomogeneous black strings 
down to D = 10 or Z) = 11, the maximum dimensions proposed in quantum 
theories of gravity, is currently unclear. 

We conclude this discussion with the caveat that the picture advocated 
here is the simplest one compatible with the numerical data currently avail- 
able. We may learn that the situation is considerably more complicated than 
we now expect, and it is certainly possible that there are even more exotic 
static black hole solutions waiting to be discovered. 
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